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Abstract — In this letter, we propose a novel adaptive reduced- 
rank strategy based on joint iterative optimization (JIO) of 
filters according to the minimization of the bit error rate (BER) 
cost function. The proposed optimization technique adjusts the 
weights of a subspace projection matrix and a reduced-rank 
filter jointly. We develop stochastic gradient (SG) algorithms for 
their adaptive implementation and introduce a novel automatic 
rank selection method based on the BER criterion. Simulation 
results for direct-sequence code-division-multiple-access (DS- 
CDMA) systems show that the proposed adaptive algorithms 
significantly outperform the existing schemes. 

Index Terms- reduced-rank techniques, adaptive algorithms, 
BER cost function, multiuser detection. 



I. Introduction 

Detecting a desired user in a DS-CDMA system requires 
processing the received signal in order to mitigate different 
types of interference. In the scenario of interest of this work, 
the receive filters are large, the interference is strong and 
time-varying, and the training data provided to the receiver 
is limited. In this context, reduced-rank signal processing has 
received significant attention in the past years, since it provides 
faster convergence speed, better tracking performance and an 
increased robustness against interference as compared to full- 
rank schemes operating with a large number of parameters. A 
number of reduced-rank techniques have been developed to de- 
sign the subspace projection matrix and the reduced-rank filter 
HI- 191 ■ Among the first schemes are the eigendecomposition- 
based (EIG) algorithms HI, JS], the multistage Wiener filter 
(MWF) investigated in 13] and Q, and the auxiliary vector 
filtering (AVF) algorithm considered in JSl. EIG, MWF and 
AVF have faster convergence speed compared to the full rank 
adaptive algorithms with a much smaller filter size, but their 
computational complexity is high. A strategy based on the 
joint and iterative optimization (JIO) of a subspace projection 
matrix and a reduced-rank filter has been reported in 0, 
Q, whereas algorithms with switching mechanisms have been 
considered in ||9l for DS-CDMA systems. 

Most of the contributions to date are either based on 
the minimization of the mean square error (MSE) and/or 
the minimum variance criteria ||T]-||9]- However, since the 
measurement of transmission reliability is the bit error rate 
(BER) not the MSE, they are not the most appropriate metric 
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from a performance viewpoint in digital communications, as 
reported in ifTOl - lfTSl . Design approaches that can minimize the 
BER have been reported in fTOl-in^ and are termed adaptive 
minimum bit error rate (MBER) techniques. The work in llT2l 
appears to be the first approach to combine a reduced-rank 
algorithm with the BER criterion. However, the scheme is a 
hybrid between an EIG or an MWF approach, and a BER 
scheme in which only the reduced-rank filter is adjusted in an 
MBER fashion. 

In this letter, we propose adaptive reduced-rank techniques 
based on a novel JIO strategy that minimizes the BER cost 
function. The proposed strategy adjusts the weights of both 
the rank-reduction matrix and the reduced-rank filter jointly in 
order to minimize the BER. By using multiple cycles over the 
recursions, we develop stochastic gradient (SG) algorithms for 
an adaptive implementation and introduce a novel automatic 
rank selection method with the BER as a metric. Simulation 
results for DS-CDMA systems show that the proposed algo- 
rithms significantly outperform existing schemes. 

II. DS-CDMA SYSTEM Model 

Let us consider the uplink of an uncoded synchronous 
binary phase-shift keying (BPSK) DS-CDMA system with K 
users, N chips per symbol and Lp propagation paths. The 
delays are multiples of the chip duration and the receiver is 
synchronized with the main path. The M-dimensional received 
vector is given by 



K 



!■(«) = ^kbk{i)Pk{^) + rikii) + n(z), 



(1) 



A;=l 



where M = N + Lp - 1, bk{i) £ {±1} is the i-th symbol for 
user k, and the amplitude of user fc is Ak- The M x 1 vector 
Pfc(i) = Cfchfe(i) is the effective signature sequence for user 
k, the M X Lp convolution matrix contains one-chip shifted 
versions of the spreading code of user k: 

I afe(l) \ 
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afc(l) 



where ak{m) e {il/ViV}, m = 1, . . . ,N. The channel vec- 
tor of user /c is hfc(j) = [hk,o{i) ■ ■ -hk^Lp-iii)]^, i7fe(«) is the 
inter-symbol interference (ISI), n(i) = [no{i) . . .nM-i{i)]'^ 
is the complex Gaussian noise vector with zero mean and 
E[n{i)n^ (i)] = cr^I, where is the noise variance, (.)^ and 
(.)^ denote transpose and Hermitian transpose, respectively. 
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III. Design of MBER Reduced-Rank Schemes 

In this section, we detail the design of reduced-rank schemes 
which minimize the BER. In a reduced-rank algorithm, an 
M X D subspace projection matrix S/j is applied to the 
received data to extract the most important information of the 
data by performing dimensionality reduction, where I < D < 
M. A D X 1 projected received vector is obtained as follows 



r(j) = Sgr(*), 



(2) 




Fig. 1. Structure of the proposed reduced-rank scheme 



where it is the input to a £> x 1 filter Wfc = [wi, W2, ■ ■ ■ , wd] ■ 
The filter output is given by Xk{i) = w^r(i) = 'wj^S^r{i). 
The estimated symbol of user k is given by bk{i) = 
sign{3fi[w^r(i)]}, where the operator 5R[.] retains the real 
part of the argument and sign{.} is the signum function. The 
probability of error for user k is given by 



Pe = Pixk <0)^ fixk)dXk 

J —oo 

' sign{b k{i)}^[xk{i)] 



Q 



(3) 



where Xk = sign{b k{i)}^[xk{i)] denotes a random variable, 
f{xk) is the single point kernel density estimate ifTOl which 
is given by 



f{ik) =- 



1 



X exp 



2'Kwj^S^SDWk 

{xk - sign{bk{tm[xkm^ 



(4) 



where p is the radius parameter of the kernel density estimate, 
Q{.) is the Gaussian error function. The parameters of S^) 
and Wfc are designed to minimize the probability of error By 
taking the gradient of Q with respect to and after further 
mathematical manipulations we obtain 



dPe 



cxp 



sign{6fe(i)} 
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,f. ^j^^uyyic;- (wf Sg SriWfc) 2 

By taking the gradient of (O with respect to SJ^ and following 
the same approach we have 



dPe 
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sign{6fc(i)} 
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IV. Proposed MBER Adaptive Algorithms 

In this section, we firstly describe the proposed scheme and 
MBER adaptive SG algorithms to adjust the weights of S£i(i) 
and w(i) based on the minimization of the BER criterion. 
Then, a method for automatically selecting the rank of the 
algorithm using the BER criterion is presented. 



A. Adaptive Estimation of Projection Matrix and Receiver 

The proposed scheme is depicted in Fig. [T] the projection 
matrix Soii) and the reduced-rank filter 'Wk{i) are jointly 
optimized according to the BER criterion. The algorithm has 
been devised to start its operation in the training (TR) mode, 
and then to switch to the decision-directed (DD) mode. The 
proposed SG algorithm is obtained by substituting the gradient 
terms ^ and (|6]l in the expressions Wfc(i + 1) = \Vk{i) — 
Hw§^ and Sr,(i + 1) = Soii)- ^^Sn■§if; CSl subject to the 
constraint of w|'^(i)S|^(i)S£i(i)wfc(i) = 1. Unlike prior JIO 
schemes, the proposed JIO-MBER algorithm employs multiple 
cycles over the recursions for Wfe and Sp. At each time instant, 
the weights of the two quantities are updated in an alternating 
way by using the following equations 



cxp 



2P 



sign{6fc(i)} 



(7) 
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2p2 



sign{6fe(z)} 



i-\J l-Kp 

X (r(z)wf (*) - S],(*)wi(z)wf 

(8) 

where [i^a and are the step-size values, the superscript j 
denotes the j-th iteration at the time instant, j = 1, . . . , J, 
and J is the maximum number of iterations. Expressions 
(|7|i and (O need initial values, w{. (0) and 8)^(0), and we 



scale the reduced-rank filter by 



at 



each iteration. The scaling has an equivalent performance to 
using a constrained optimization with Lagrange multipliers 
although it is computationally simpler The updated filters 
for the next time instant are given by w^(i + 1) ^ Wjf(j) 
and S^(i + 1) ^ S^(i). The proposed adaptive JIO-MBER 
algorithm is summarized in table J] 

Since the BER is employed as a design criterion, there is no 
guarantee that the algorithms proposed will obtain the global 
minimum solution IIO]. The proposed JIO-MBER technique 
and the other existing BER-driven algorithms can only have 
their convergence guaranteed to a local minimum. In |[T6l . 
the work provides a general proof for the convergence of 
alternating optimization algorithm. By extending the work in 
|[T6l . we can obtain local minimum solutions for the projection 
matrix and the reduced-rank filter of the proposed JIO-MBER 
algorithms. 
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TABLE I 

Proposed adaptive JIO-MBER algorithms 



Initialize w^(0) and S^(0). 

Set step-size values ^-w and liSo "'^ iterations J. 

for each time instant i do 
for j from 1 to J do 

Compute w^"*" (i) and (i) using (7) and (8). 

Scale the using < , 

end cycles 

After the J cycles, obtain w^(i + 1) <— w^(i) 

and S|^(^ + 1) <— Sj^(i) for the next time instant. 



B. Computational Complexity of Algorithms and Delay Issues 

We describe the computational complexity of the proposed 
JIO-MBER adaptive algorithm in DS-CDMA systems. In 
Table HIl we compute the number of additions and multi- 
plications to compare the complexity of the proposed JIO- 
MBER algorithm with the conventional adaptive reduced-rank 
algorithms, the adaptive least mean squares (LMS) full-rank 
algorithm based on the MSE criterion ifTSi and the SG full- 
rank algorithm based on the BER criterion ifTol. Note that 
the MWF-MBER algorithm corresponds to the use of the 
procedure in ||3] to construct Soii] and ^ with J = 1 to 
compute w[z]. In particular, for a configuration with N = 31, 
D = 6 and Lp ~ 3, the number of multiplications for 
the MWF-MBER and the proposed JIO-MBER algorithms 
are 8377 and 1262, respectively. The number of additions 
for them are 5920 and 962, respectively. Compared to the 
MWF-MBER algorithm, the JIO-MBER algorithm reduces the 
computational complexity significantly. 

TABLE II 

Computational complexity of Algorithms. 



C. Automatic Rank Selection 

The performance of reduced-rank algorithms depends on the 
rank D, which motivates automatic rank selection schemes to 
choose the best rank at each time instant [3], Q, Q. Unlike 
prior methods for rank selection, we develop a rank selection 
algorithm based on the probability of error, which is given by 



PDit) = Q 



sign{6fc(i)}5i[xf(i)] 



(9) 



where the receiver is subject to wj^ S^Sn'Wk = 1. 
For each time instant, we adapt a reduced-rank filter 
Wfc(j) and a projection matrix S£)(i) with the max- 
imum allowed rank Dmax, which can be expressed 



as wfc(i) 
where Drr 



SM,l(«) 

and i?n 



are the 



Xi),...,WD„,,,^{i)] and 

... Sl,A„ax(«) 



minimum 



and maximum 



ranks allowed for the reduced-rank filter, respectively. For 
each symbol, we test the value of rank D within the range. 



namely, 



< D < D„ 



substitute the filter 



matrix S^(i) 



.Wfc(i) = 

Sl,l(«) 



(. For each tested rank, we 

[uii(i), . . . , ^^(i)]^ and the 
... si,_d(«) 

into ^ to 

sm,i(0 • ■ • SM,D{i) 
obtain the probabiUty of error Poii)- The optimum rank can 
be selected as 



-Dopt(i) = arg min 

DG{r'™,„,...,-D„ 



(10) 



Algorithm 


Number of operations per symbol 


Multiplications 


Additions 


tull-Rank-LMS 


'IM + 1 


'2M 


Full-Rank-MBER 


4Af + l 


AM - 1 


MWF-LMS g) 


DM'^ - 


DM^ - 




+2DM + 4D + 1 


+3D-2 


EIG dl 






JIO-LMS (D 


SDM + M 


2DM + M 


+3D + 6 


+4D - 2 


MWF-MBER QU 


(D + 1)M^ 


(D - l)A/2 
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+MLp + 1 


JIO-MBER 


6MDJ + 5DJ 


5MDJ + DJ 
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-MJ - J 



Note that generally the delay associated with receive pro- 
cessing of filters is quite small as compared to the delay of 
channel decoding used in some systems such as maximum a 
posteriori (MAP), sum-product algorithm (SPA) and others for 
turbo and LDPC codes. We can find that the delay introduced 
by the receive filter is much lower than that of the channel 
decoder in [[TtI- We can see that the decoder has a much 
bigger impact on the delay and actually dominates the delay 
in the receiver ifTSl . For example, in the case of 64 kbit/s 
UMTS transmission, the encoded broadcast service requires 
more than 5 s system delay, and the uncoded systems such as 
conferencing and telephony require 150 ms delay. 



Note that after we obtain the projection matrix with rank 
Dmax, we generate the M x D projection matrix by grouping 
the columns which are from the first column to the D-th 
column. The proposed auto rank selection mechanism only 
requires a modest increase in the complexity of the JIO- 
MBER algorithm with a fixed rank. The number of multi- 
plications and additions for the JIO-MBER algorithm with 
the auto rank mechanism are (6A/ + 5)Dmax + A/ + 11 and 
(5A/ + l)Dmax — M — 1, respectively. In addition, a simple 
search over the values of x^{i) and the selection of the terms 
corresponding to Dopt and Pn^^t (*) ^^e performed. 



V. Simulations 

In this section, we evaluate the performance of the pro- 
posed JIO-MBER reduced-rank algorithms and compare them 
with existing full-rank and reduced-rank algorithms. Monte- 
carlo simulations are conducted to verify the effectiveness 
of the JIO-MBER adaptive reduced-rank SG algorithms. The 
DS-CDMA system employs Gold sequences as the spread- 
ing codes, and the spreading gain is iV = 31. The se- 
quence of channel coefficients for each path is hkj{i) = 
Pk.fCtkj{i){f — 0,1,2). All channels have a profile with 
three paths whose powers are pk.o = dB, pk.i = —7 
dB and pfc 2 = —10 dB, respectively, where akj{i) is 
computed according to the Jakes model. We optimized the 
parameters for the proposed and conventional algorithms based 
on simulations. In this work, we focus on the case with fixed 
step-size values and we note that, motivated by the work in [?], 
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low-complexity variable step-size mechanisms with the error 
probability as a metric can be further developed. The initial 
full rank and reduced-rank filters are all zero vectors. The ini- 
tial projection matrix is given by 8)5(0) = [I_d, 0£)x(M-n)]"^- 
The algorithms process 250 symbols in TR and 1500 symbols 
in DD. We set p = 2a. 



N=31 .SNR=15 dB, K-5, 




500 1000 1500 

number ot received symbcis 

Fig. 2. BER perforaiance versus number of received symbols for tlie JIO- 
MBER reduced-ranlc algoritiims and tiie conventional schemes. JIO-MBER 
(D=8, J=5): ^t„, = = 0.005. JIO-MBER (D auto, J=5): ^u, = 113 ^ = 
0.005. JIO-MBER (D=8, J=l): fiyj = ^J.SD = 0-105. JIO-MBER (D auto, 
J=l): = = 0.16. EIG-MBER (D=8): ^tu. = 0.215. LMS full- 

rank: /i„ = 0.105. MBER full-rank: /i» = 0.05. MWF-MBER (D=8): 

At™ = 0.05. (Dmin = 3, Dmax = 20) 

Fig. |2] shows the bit error rate (BER) performance of 
the desired user versus the number of received symbols for 
the JIO-MBER adaptive SG algorithms and the conventional 
schemes. We set the rank D = 8, K = 5, SNR = 15dB and 
fj^Ts = 5 X 10"^ We can see that the JIO-MBER reduced- 
rank algorithms converge much faster than the conventional 
full rank and reduced-rank algorithms. For the group of JIO- 
MBER adaptive algorithms, the auto-rank selection algorithms 
outperform the fixed rank algorithms. Fig. |3] illustrates the 
BER performance of the desired user versus SNR and number 
of users K. In particular, the JIO-MBER algorithm using 
D = 8 with ,7 = 1 iteration can save up to over 6dB and 
support up to six more users in comparison with the MWF- 
MBER algorithm using £) = 8, at the BER level of 2 x 10"^ 

VI. Conclusions 

In this paper, we have proposed a novel adaptive MBER 
reduced-rank scheme based on joint iterative optimization of 
filters for DS-CDMA systems. We have developed SG-based 
algorithms for the adaptive estimation of the reduced-rank 
filter and the projection matrix, and proposed an automatic 
rank selection scheme using the BER as a criterion. The 
simulation results have shown that the proposed JIO-MBER 
adaptive reduced-rank algorithms significantly outperform the 
existing full-rank and reduced-rank algorithms at a low cost. 
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